Research in Astron. Astrophys. Vol.0 (201x) No.O, 000-000 
http:/ /www. raa-journal. org http:// www. iop. org/ journals /raa 



Research in 
Astronomy and 
Astrophysics 



Bianchi Type V Universe with Bulk Viscous Matter and Time 
Varying Gravitational and Cosmological Constants 



Prashant S. Baghel* and J. P. SingbJ 

Department of Mathematical Sciences, A. P. S. University, Rcwa 
* ps.baghel@yahoo.com, ^ jpsinghmp@lycos.com 



486003, India; 



o 

(N 



(N 



> 
in 

(N 

in 

in 
o 

(N 



X 



Abstract Spatially homogeneous and anisotropic Bianchi type V space-time with bulk 
viscous fluid source and time varying gravitational constant G and cosmological term A 
are considered. Coefficient of bulk viscosity £ is assumed as a simple linear function of 
Hubble parameter H (i.e. ( = + CiH, where Co and £i are constants). The Einstein field 
equations are solved explicitly by using a law of variation for the Hubble parameter, which 
yields a constant value of deceleration parameter. Physical and kinematical parameters of the 
models arc discussed. The models arc found to be compatible with the results of astronomical 
observations. 

Key words: Bulk viscosity; Bianchi V; Constant deceleration parameter; Variable cosmo- 
logical term. 



1 INTRODUCTION 

In Einstein theory of Gravity Newtonian gravitational constant G and Cosmological term A are con- 
sidered as fundamental constants. The Newtonian constant of gravitation G plays the role of coupling 
constant between geometry of space and matter in Einstein field equ ations . In an evolving universe it 
is natural to look at this constant as a function of time. Ever since iDirad (|l937h first considered the 
pos sibility of varia ble G, there has been numerous modifications of general relativity to allow a variable 
G ( Wesson! . Il980h. Neverthel e ss, these theories have not gained wid e acceptanc e. However, recently a 
modification (jAbdel-Rehmanl . Il992t iBeeshaml Il986at iBermanl . Il99lt ILaul . Il985h has been proposed in 
Ein stein field equations t r eating G and A as coupling variables within the framework of general relativ- 
ity. ICanuto and Narlikarl (|l980f) have s hown that the G varyi n g cosmology is consistent with whatever 
cosmological observations are available. lAbdel-Rehnianl (|l990l) ; feeeshanl (ll986al ); iLevitl (|l980h have dis- 
cussed the possibility of an increasing G . 

Cosmological models with a cosmological constant are currently serious candidates to describe the ex- 
pansion history of the universe. The huge difference between the small cosmological constant inferred from 
observations and vacuum energy density re s ulting from quantum field theories has been for along time, a 
diffic ult and fascinating problem ( Carneircl |2003[ [kraus sl ll995HSahni and Starobinskvl. 120001; IWeinberd . 
19891) for cosmologists and field theo r y researchers. A wide range of observations ( Garnavich et all 1998 ; 
Perlmutter et al. 119991: iRiess et all . l2004t ISchmidtl . Il998h suggest that universe possesses a non-zero 
cosmological constant. At the same time, theories are increasingly exploring the possibility that this pa- 
rameter is a dynamical one whose effective value in the early universe may have been quite different from 
the one we measure today. One possible explanation for a small A term is to assume that it is dynamically 
evolving and not constant, that is, as the universe evolves from an earlie r hotter and denser epoch, the 
effective cosmological term also evolves and decre a ses to its p r esent value ( Coble etail Il997t iFrees et ail 
1987t iFrieman etafl Il995t lOzer and Tahal . Il987l ). iBertolamil (|l986l ) obtained the time-dependent G and 
A solutions. Cosmological models with variable G and A term have been studied by a number of authors 
for homogeneous isotropic (lAbdel-Rehmanl . Il990b lAbdussattar and Vishwakarma . 1997 ; Arbabl. I l998a ; 



Beesham L 1986bt Kalligas et al. , 1992 ) and an i sotrop ic ( Bali and Tinker L 20091; Pradhan et all 120051 ; iSaha 
20051: ISineh and BeeshamL 1201ft ISingh et all 120081 lYadavl l2010l lYadav et alll2012h space-times 



The investigations of relativistic cosmological models usually has the cosmic fluid as perfect fluid. 
However, these models do not incorporate dissipative mechanisms responsible for smoothing out ini- 
tial anisotr opics. I t is b elieved that during neutrinos dec oupling , the matter behaved like a vis- 
cous fluid (iKlimekl. [ 19761 ) in the early stages of evolution. IColevI (|l990t ) studied Bianchi V viscous 
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fluid cosmologi cal models for barotropic fluid distribution. iHeller and Klimekl (| 19751) : iMisnerl (ll967f ) 



iMurphvl (11973) have investigated t he role of viscosity in avoiding the initial big bang singular 



ity. Padmanabhan and ChitreT ( 19871 ) investigated the effect of bulk viscosity on the evolution of the 



univer se at large. They have shown that bulk viscosity lea ds to the infla ti onary like solutions. ISahal 
( 2005[ ) discussed Bianchi type I universe with viscous fluid. iPradhan et al.l ( 2004 ) investigated viscous 



fluid cosmological models in LRS Bianchi type V universe with var ying A. Bulk v i scous cosmologi- 
cal m o dels with time-d e pende nt G and A term have been studied bv lArbabl ( 1998bl ): iBali and Tinker! 



cai models witn time-de pende nt U ana a term nave been studied bv lArba b ^iyy»bjj; man and unitei 
(|2009[ ): iBeesham etal] (|2000h. where bulk viscosity is taken as a power function of energy density. 



Recently Singh and BaghelT ( 2010l) investigated Bianchi type V cosmological models with bulk viscos 



ity, where the coefficient of bulk viscosity is assumed to be a power function of energy density p or 
volume expansion 9. 

Among the physical quantities of interest in cosmology the deceleration parameter q is currently a 
serious candidate to describe the dynamics of the universe. The prediction of standard cosmology, that 
the universe at present is decel erating is contradictory to the recent observational evidences of the high 



the universe at present is decelerating is con tradictory to the recent obs ervational evidences ot the high 
redshift of type la supernovae ( Knori 120031 : iRiess etal 1. 119981: iTornV et al. 1. 120031) . Observations reveal 



that instead of slowing down, the expanding universe is speeding up. Models with constant decelera- 
tion parameter have received considerable att ention recently. The law of variation for Hubble param- 
eter was initially proposed by iBerma 3 (119831 ) for FRW models that yields a constant value of decel- 



eration parameter. Cosmologi cal models with such a law of variation for Hubble parameter have bee n 
studied by number of authors (iMaharai and Naidool . Il993t iReddv et al. |, l2007t ISingh and Kumarl . I2009D . 
Recently ISingh and" Baghell (|2009aU b() proposed a similar law for the Hubble parameter and generated the 
solution for Bianchi type V space-time in general relativity. According to the proposed law the relation 
between the Hubble parameter H and average scale factor R is given by 

H = /3iT m , (1) 

where /3 > and m > are constants. For this variation law the deceleration parameter q comes out to 
be constant, i.e. 

q = m-l. (2) 

For m > 1, the model represents a decelerating universe and m < 1 corresponds to the accelerating phase 
of the universe. When m = 1, we obtain H ~ i a nd q = 0. Therefore galaxies move with constant speed 
and model represents anisotropic Milne universe (Landsbe rg and Evansl . fl977l ) for m = 1. For m = 0, we 



get H = (3 and q = — 1. Thus the observed Hubble parameter is true constant equal to its present value 
Hq and model represents accelerating phase of the universe. 

The releva nce of the study of Bianch i type V cosmological models has already been discussed in our 
earlier papers (ISingh and Baghel l2009al ibh. where we studied the viscous fluid models in some details. 



As a natural sequel to that study here we incorporate time varying G and A term in the bulk viscous 
Bianchi type V models and solve the coupled field equations exactly. In thi s paper we have ta ken the 
coefficient of bulk viscosity as simple linear function of Hubble parameter H ( Meng et al.l . 120071 ) . 



2 METRIC AND FIELD EQUATIONS 

We consider Bianchi type V space-time in orthogonal form represented by the line-element 

ds 2 = -dt 2 + A 2 (t)dx 2 + e 2ax {B 2 {t)dy 2 + C 2 (t)dz 2 } , (3) 

where a is a constant. We assume that the cosmic matter is represented by the energy- momentum tensor 
of an imperfect bulk viscous fluid 

T ij = (P + P)viVj + pg l0 , (4) 
where p is the effective pressure given by 



satisfying linear equation of state 



P=P-(v i ;i, (5) 



V = up. (6) 



Here p is the equilibrium pressure, p is the energy density of matter, £ is the coefficient of bulk viscosity and 
v l is the flow vector of the fluid satisfying v^v 1 = — 1. The semicolon stands for covariant differentiation. 
On thermodynamical grounds bulk viscous coefficient £ is positive, assuring that the viscosity pushes 
the dissipative pressure p towards negative values. But correction to the thermodynamical pressure p 
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due to bulk viscous pressure is very small. Therefore, the dynamics of cosmic evolution does not change 
fundamentally by the inclusion of viscous term in the energy-momentum tensor. 

The Einstein fields equations with time-dependent cosmological term A and gravitational constant G 

are 

Rij - = -^G{t)T l3 + A(t) 9ij . (7) 

For the metric ^ and matter distribution Q in comoving system of coordinates (vi ~ —fif), the field 
equations (|Tj) yield 

o ^ . " 2 C A CA 

S * Gp - A =Al-C-A~CA> (9) 

^--A=^-4-|-4i, do) 

A 2 A B AB v ' 
o n 3 « 2 AB BC AC 

8nGp + A = -^ + AB + BC + AC> (U) 
n 2A B C 

° = ^-B-C> (12) 

where an overhead dot (.) denotes ordinary differentiation with respect to cosmic time t. Due to contracted 
Bianchi identity, divergence of Einstein tensor = R%j — \R\gij is zero and we get 

SttG jp+( P + 3 5) ^ + | + gj|+87rG + A = 0. (13) 

The usual conservation equation T( = splits above equation into 

P+(p + p)(j + § + ^\=0 (14) 

and 

87rpG + A = 0. (15) 

From p5p. one concludes that, when A is constant or zero, G turns out to be constant for non-zero energy 
density. 

We define the average scale factor R for Bianchi V space-time as R 3 = ABC. From equations (|5])- (fTU)) 
and ([12"]) , we obtain 

A R 

A = R> ™ 

(17) 

(18) 

k being constant of integration. On integration, equations (llGD - (jl8p give 

A = miR, (19) 

B = m 2 Rexp (-k J - J J , (20) 

C = m 3 Rexp(kJ , (21) 
where mi, mi and 7713 are constants of integration satisfying m\m-2jriz = 1. 
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In analogy with FRW universe, we define a generalized Hubble parameter H and generalized decel- 
eration parameter q as 

H = ^ = ^(H 1 + H 2 + H 3 ) (22) 

and 

9 = -l-ff (23) 

where Hi = 4, -H2 = 5 an d = § are directional Hubble's factors along x, y and z directions 
respectively. 

We introduce volume expansion 9 and shear a, as usual 

9 = u?j and er 2 = -oyer 13 , 
fj lJ being shear tensor. For the Bianchi V metric, expressions for 9 and a come out to be 

(24) 



3R 
~R 



and 



(25) 



k 

We can express equations (|5))- ([TTj) and (fTl)) in the terms of R, H, q and a as 

a 2 

8irGp-A = (2q-l)H 2 -(T 2 + —, (26) 

R z 

3a 2 

8ttG P + A = 3H 2 -a 2 (27) 

p + 3(p + p)iJ = 0. (28) 

It is to note that energy density of the universe is a positive quantity. It is believed that at the early 
stages of evolution, when the average scale factor R was close to zero, the energy density of the universe 
was infinitely large. On the other hand, with expansion of the universe i.e. with increase of R , the energy 
density decreases and an infinitely large R corresponds to p close to zero. In that case from (|27[) . we 
obtain ^ — > 1, where p v = and p c = For A > 0, p < p c . Also from (f2"T)) . we observe that 

< §2 < 5 arid < 8 ^ p < | for A > 0. Thus a positive A restricts the upper limit of anisotropy 
whereas a negative A will increase the anisotropy. 
From (gU) and (gTJ), we get 

^ = A + 12ttG(9 - AttG{ P + 3p) - 2a 2 - °— , (29) 
etc o 

which is the Raychaudhuri equation for the given distribution. We observe that for negative A and 
the absenc e of viscosity, the univers e will always be in decelerating phase provided the strong energy 
conditions ([Hawking and Ellisl . Il975h hold whereas in the presence of viscosity, positive A will slow down 



the rate of decrease of volume expansion. Also, a = —o9 implying that a decreases in an evolving universe 
and for infinitely large value of i?, a becomes negligible. 

3 SOLUTION OF THE FIELD EQUATIONS 

Integrating ([T]), we obtain 

R = (mf3t + h)™ for m ^ (30) 

and 

R = cxp{(3{t - t )} for m = 0, (31) 
where t\ and to are constants of integration. From (|T9"]) - (|21[) with the use of ((30]) . we obtain 



A = mi(m/3t + ti)- 
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B = rri2 (m/3t + 1 \ ) m exp 



C = m 3 (mf3t + ti) m exp ■ 



-k{m/3t + ti 



P(m - 3) 
k(mj3t + ti)- 



/?(m - 3) 

For this solution, the metric ([3]) assumes the following form after suitable transformation of coordinates 

ds 2 = -dT 2 + (m/3T)^dX 2 + [mpT)% exp LaX - T) ~™~ ) dY 2 

P(m - 3) J 

+(m/3T)* exp ^2aX + } (32) 

Equations (JT9J) - (J2T|) together with dSTJ give 

A = miexp{/3(i - t )}, 

i? = m 2 cxp|/3(t-to) + ^e- 3 ^-^ 

C = m 3 exp U(t - t Q ) - ^~ mt ~ to) 
The line-element ([3]) for this solution can be written as 

+ exp (^2aX + 2/3T - i^e" 3 ^ rfZ 2 . (33) 

4 DISCUSSION 

To determine the coefficient of bulk viscosity, £ is assumed to be a simple linear function of Hubble 
parameter H ( Meng et all l2007h . i.e. 

C = Co + Citf, (34) 
where Co (> 0) and Ci are constants. For this choice, equation (j28|) reduces to 

p + 3(l + Lj)H P = 9((o + (iH)H 2 . (35) 

We discuss the models for m ^ and m = 0. 

4.1 Cosmology for m 7^ 

For the model (|3"2"j) , average scale factor i? is given by 

fl=(m/3T)™. (36) 
Volume expansion 0, Hubble parameter H and shear scalar a for this model are 

= 3H=JL, (37) 

a = k(m/3T)^. (38) 

Using fl37j) in (J34J) and ([35]), we obtain 

C = Co + ^, (39) 

= 9Co 9Ci a 

P m(3 + 3w - m)T (3 + 3w - 2ra)m 2 T 2 ^3(1+^)/™ ' ^ > 
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Time t 



Fig. 1 Variation of matter energy density p with cosmic time t in case of m ^ 0. 



where a is an integration constant. Gravitational constant G and cosmological term A are obtained as 

mT 2 {a 2 (m^T)£ + fc 2 } - (mf3T)% 



G= — 



1 

4tt 



2m + 3(l-u J ) 



(41) 



A = 18 



Ci 



o 



Ci 



9T 



31 ! I j» 



(3 + 3oj — m)mT (3 + 3cj — 2m)m 2 T 2 
(m/3T)™ - mT 2 |a 2 (m/3T)^ + fc 2 } 



(^fe) ^ + (as^zj) Ti - ma(l + ^ 
3 



+ - 



T 
A- 2 



2m— 3a; — 3 

3a 2 



2^2 



(mf3T)- (m(3T)- 



(42) 



We observe that the model is not tenable for to = 3. The model has singularity at T = 0. At T = 0, 
p, A, C, 0, <7 all diverge whereas G becomes a constant provided m > 3. In the limit of large times, 
p, 9, a become zero, C = Co an d G — > oo, A — ^constant for m > 3. Also for T — > oo, er/# — > when 
to < 3. Therefore the model approa c hes iso tropy asymptotically. For large values of T the model becomes 
conformally flat ( Singh and Baghel . l2009bT) . The integral 



dt 



R(t) /3(m - 1) 



(mf3T) 2 



is finite provided ra/1. Therefore particle horizon exists in the model. It is to note that for to 
model does not have horizon. 



(43) 
1, the 



4.2 Cosmology for m = 

We now discuss the model (|3"3"]l . Average scale factor R, expansion scalar 9, Hubble parameter H shear 
a and deceleration parameter q are given by 

R = e? T , (44) 

9 = 3H = 3/3, (45) 
a = ke~^ T , (46) 
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Fig. 2 Variation of matter energy density p with cosmic time t in case of rn = 0. 



q = -\. (47) 
We obtain bulk viscous coefficient £ and matter density p as 

C = Co + Ci/3, (48) 

p = 3/3(Co + Ci/3) + &e- 3W +") T , (49) 
where 6 is a constant of integration. Expressions for gravitational constant G and cosmological term A 

(fc 2 + a 2 e^ T ) 



are 



G = — 



47rfe(l 



=3(1- 



A = 3/3 + 





* 2 - 




a* 




e 6/3T 


V 1 + w J 


e 2f3T 



6/3(Co 



Ci/3)(fc 2 



a e 



4/3T\ 



(50) 



(51) 



&(l + w)2 e 3(i-")/3T 

The model has no initial singularity. The expansion scalar 6 is constant throughout the evolution of the 
universe and therefore the model represents uniform expansion. Also q = — 1 shows that the expansion of 
the model is accelerating with constant speed rate. At T = 0, R, a, £, p. G and A all are constant. The 
energy density decreases as time increases and becomes constant at late times. Also G and A become 
infinite for very large values of T whereas they tend to zero for to < — | at late times. We observe that 
gravitational constant G to be positive for b < wherea s G is negat i ve for b > 0. The possibility of 
negative gravitational constant G has been discussed bv IStarobinskvl ( 198lh concluding that effective 
gravitational constant may have changed sign in the early universe. As T — > oo, the ratio a/6 becomes 
zero. Therefore the model approaches isotropy for large values of T. 



5 CONCLUSION 



In this paper, we have investigated spatially homogeneous and anisotropic Bianchi type V space-time with 
bulk viscous matter and time-dependent gravitational constant G and cosmological term A in general 
relativity. The field equations have been solved exactly by using a law of variation for the generalized 
Hubble parameter H. Two universe models have been obtained and physical behaviour of the models is 
discussed. Coefficient of bulk viscosity £ is assumed to be a simple linear function of Hubble parameter H 
(i.e. £ = ( + (iH, where £ and Ci are constants). For ( = = d we recover perfect fluid models. In case 
of cosmology for m^O, the universe starts from a singular state whereas cosmology for m = follows a 
non-singular start. We observe that the presence of bulk viscosity increases the value of matter density. 
For the models obtained a/6 — > as T - >oo. Thus the models approach isotropy at late times. From 
([2]) one concludes that for m > 1, the model represents a decelerating universe whereas for < m < 1, it 
gives rise to an accelerating universe. When m = 1, we obtain H = i and q = 0, so that every galaxy 
moves with constant speed. Behaviour of the matter energy density p and vacuum energy density A with 
respect to cosmic time are plotted in figures [U [21 El Ql 
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Fig. 3 Variation of vacuum energy density A with cosmic time t in case of m 
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Fig. 4 Variation of vacuum energy density A with cosmic time t in case of m 



References 

Abdel-Rehman, A.M.M. 1990, Gen. Relativ. Gravit., 22, 655 

Abdel-Rehman, A.M.M. 1992, Phys. Rev. D, 45, 3497 

Abdussattar and Vishwakarma, R.G. 1997, Class. Quant. Grav., 14, 945 

Arbab, A.I. 1998a, Gen. Relativ. Gravit., 29, 391 

Arbab, A.I. 1998b, Gen. Relativ. Gravit., 30, 1401 

Bali, R. and Tinker, S. 2009, Chin. Phys. Lett., 26, 029802 

Beesham, A. 1986a, Int. J. Theor. Phys., 25, 1295 

Beesham, A. 1986b, Nuovo Cimento B, 96, 17 

Beesham, A., Ghosh, S.G. and Lombard, R.G. 2000, Gen. Relativ. Gravit. ,32,471 
Berman, M.S. 1983, Nuovo Cim. B, 74,182 
Berman, M.S. 1991, Gen. Relativ. Gravit., 23, 465 
Bertolami, O. 1986, Nuovo Cim. B, 93, 36 



Bianchi type V universe with bulk viscous matter and time varying constants. 



9 



Canuto, V.M. and Narlikar, J.V. 1980, Astrophys. J., 6, 236 
Carneiro, S. 2003, Int. J. Mod. Phys. D, 12, 1669 

Coble, K., Dodelson, S. and Frieman, J.A. 1997, Phys. Rev. D, 55, 1851 

Coley, A. A. 1990, Gen. Relativ. Gravit., 22, 3 

Dirac, P.A.M. 1937, Nature, 139, 323 

Frees, K. et al. 1987, Nucl. Phys. B., 287, 797 

Frieman, J.A. et al. 1995, Phys. Rev. Lett., 75, 2077 

Garnavich, P.M. et al. 1998, Astrophys. J., 509, 74 

Hawking, S.W. and Ellis, G.F.R. 1975, The Large Scale Structure of Space-time, Cambridge University 
Press, p. 88 

Heller, M. and Klimck, Z. 1975, Astrophys. Space Sci., 33, 37 

Kalligas, O., Wesson, PS. and Everitt, C.W.F. 1992, Gen. Relativ. Gravit., 24, 351 
Klimek, Z. 1976, Nuovo Cim. B, 35, 249 
Knop, R.A., et al. 2003, Astrophys. J., 598, 102 
Krauss, L.M. et al. 1995, Gen. Relativ. Gravit., 27, 1137 

Landsberg, P.T. and Evans, D.A. 1977, Mathematical Cosmology, Clarendon, Oxford, p. 101 
Lau, Y.K. 1985, Aust. J. Phys., 38, 547 
Levit, L.S. 1980, Lett. Nuovo Cimcnto, 29, 23 

Maharaj, S.D. and Naidoo, R. 1993, Astrophys. Space Sci., 208, 261 

Meng, X.H., Ren, J. and Hu, M.G. 2007, Commun. Thcor. Phys., 47, 379 

Misner, C.W. 1967, Nature, 214, 40 

Murphy, J.L. 1973, Phys. Rev. D, 8, 4231 

Ozer, M. and Taha, M.O. 1987, Nucl. Phys. B., 287, 776 

Padmanabhan, T. and Chitre, S.M. 1987, Phys. Lett. A, 120, 433 

Perlmutter, S. et al. 1999, Astrophys. J., 517, 565 

Pradhan, A., Yadav, L. and Yadav, A.K., 2004, Czech. J. Phys., 54, 487 

Pradhan, A., Yadav, A.K. and Yadav, L., 2005, Czech. J. Phys., 55, 503 

Reddy, D.R.K., Naidoo, R.L. and Adhav, K.S. 2007, Astrophys. Space Sci. 307, 365 

Riess, A.G. et al. 1998, Astron. J., 116, 1009 

Riess, A.G. et al. 2004, Astrophys. J., 607, 665 

Saha, B. 2005, Mod. Phys. Lett. A, 20, 2127 

Saha, B. 2006, Astrophys. Space Sci., 302, 83 

Sahni, V. and Starobinsky, A.A. 2000, Int. J. Mod. Phys. D, 9, 373 
Schmidt, B.P. 1998, Astrophys. J., 507, 46 

Singh, CP. and Beesham, A. 2010, Int. J. Mod. Phys. A, 25, 3825 

Singh, CP. and Kumar, S. 2009, Int. J. Theor. Phys., 48, 2401 

Singh, J. P. and Baghel, PS. 2009a, Electronic J. Theor. Phys., 6, 85 

Singh, J.P. and Baghel, PS. 2009b, Int. J. Theor. Phys., 48, 449 

Singh, J.P. and Baghel, PS. 2010, Int. J. Thcor. Phys., 49, 2734 

Singh, J. P., Pradhan, A. and Singh, A.K. 2008, Astrophys. Space Sci., 314, 83 

Starobinsky, A.A. 1981, Sov. Astron. Lett., 7(1), 36 

Tonry, J.L., et al. 2003, Astrophys. J., 594, 1 

Weinberg, S. 1989, Rev. Mod. Phys., 61, 1 

Wesson, P.S. 1980, Gravity, Particles and Astrophysics (Reidel, Dordrecht). 
Yadav, A.K., 2010, Int. J. Theor. Phys., 49, 1140 

Yadav, A.K., Pradhan, A. and Singh, A.K. 2012, Astrophys. Space Sci., 337, 379 



